Formation Rates of Black Hole Accretion Disk Gamma-Ray Bursts

While many models have been proposed for GRBs, those currently favored are all based upon the formation of and/or rapid accretion into stellar mass black holes. We present population synthesis calculations of these models using a Monte Carlo approach in which the many uncertain parameters intrinsic to such calculations are varied. We estimate the event rate for each class of model as well as the propagation distance for those having significant delay between formation and burst production, i.e., double neutron star (DNS) mergers and black hole-neutron star (BH/NS) mergers. For reasonable assumptions regarding the many uncertainties in population synthesis, we calculate a daily event rate in the universe for i) merging neutron stars: ~100/day; ii) neutron-star black hole mergers: ~450/day; iii) collapsars: ~10,000/day; iv) helium star black hole mergers: ~1000/day; and v) white dwarf black hole mergers: ~20/day. The range of uncertainty in these numbers however, is very large, typically two to three orders of magnitude. These rates must additionally be multiplied by any relevant beaming factor and sampling fraction (if the entire universal set of models is not being observed). Depending upon the mass of the host galaxy, half of the DNS and BH/NS mergers will happen within 60kpc (for a Milky-Way massed galaxy) to 5Mpc (for a galaxy with negligible mass) from the galactic center. Because of the delay time, neutron star and black hole mergers will happen at a redshift 0.5 to 0.8 times that of the other classes of models. Information is still lacking regarding the hosts of short hard bursts, but we suggest that they are due to DNS and BH/NS mergers and thus will ultimately be determined to lie outside of galaxies and at a closer mean distance than long complex bursts (which we attribute to collapsars).Comment: 57 pages total, 23 figures, submitted by Ap

A line-binned treatment of opacities for the spectra and light curves from neutron star mergers

The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion-opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and that eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion-opacity formalism. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semi-relativistic approaches.Comment: 27 pages, 22 figures. arXiv admin note: text overlap with arXiv:1702.0299